Quantum buses and quantum computer architecture based on quantum dots 



in 

O 

o 

> 
O 

oo 



Irene D'AmiccQ 

Department of Physics, University of York, York YOlO, United Kingdom 
(Dated: February 2, 2008) 

We propose a quantum computer architecture based on quantum dots both for short distance and for long 
distance communication/computation. Our scheme exploits the natural characteristics of self-assembled quan- 
tum dots and it is scalable. It is centered on the idea of a quantum bus based on semiconductor self-assembled 
quantum dots. This allows for transmission of qubits between the different quantum registers, and could be 
integrated in most of the present proposal for semiconductor quantum dot-based quantum computation. Our 
proposal exploits the peculiar properties of relatively short spin-chains, and advantages and disadvantages of 
two possible implementations, both based on spin-chain global dynamics, are discussed in details. A clear ad- 
vantage of the scheme is to avoid the use of microcavities for long distance communication between different 
elements of the quantum computer. In this respect our scheme is comparatively faster than hybrid quantum 
dot-microcavity schemes. 
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I. INTRODUCTION 

In the past years, self-assembled semiconductor quantum 
dotsi have been considered one of the most promising hard- 
ware for envisaging quan timi inforrriat ion/computation (QIC) 
schemes (see e.g. Ref. 2345678). The key feature of 
this hardware is the possibility of optically drive quantum cal- 
culations by multicolor trains of laser pulses, overcoming in 
this way the relatively fast decoherence times typical of solid 
state systems. Another important characteristic is the fact that, 
since self-assembled quantum dots (QDs) tend to be naturally 
produced in large arrays, such a hardware could in principle 
solve the problem of scalability, which affects other types of 
implementations . 

Many QD-based QIC schemes consider as quantum arrays 
ensembles of vertically stacked quantum dots^ -'-"^'^'^, which 
spontaneously form in multilayer QD structures'. However 
following the initial enthusiasm, it has become clear that this 
choice solves only partially the scalability problem. Experi- 
ments show in fact that the maximum length of such vertical 
arrays is of the order of ten quantum dots only. Clearly this 
rises the problem of how to perform complex computations or 
equivalently, how to communicate between different aiTays. A 
possible solution for communicating between far away quan- 
tum dots/arrays is to embed the system in a microcavity and 
then to use the cavity mode to transmit qubits of information 
between well separated quantum dots. This exploits the po- 
tential conversion of trapped excitons into flying qubits and 
vice-versa^ '^. Among the problematic issues related to this 
scheme are (i) to experimentally embed self-assembled arrays 
of quantum dot in a high Q microcavity, (ii) to carefully posi- 
tion the dots away from the cavity mode nodes in order to pre- 
serve QD-cavity coupling, (iii) to engineer the structure such 
that it is still possible to address the QDs by laser fields and 
possibly couple each QD with a static external electric fields, 
as required, e.g., in Ref. I2l3l7i Experimentally this is very 
demanding. 

To lift these complex requirements, we propose to commu- 
nicate between far away vertical arrays of QDs by using quan- 
tum buses made by in-plane chains of quantum dots. To im- 



plement this, our proposal exploits the structural differences 
and peculiar characteristics of in-plane versus vertical arrays 
of self-assembled QDs. 

The first part of the paper describes a possible global ar- 
chitecture for QD-based quantum computers. In the second 
part we will analyze different possibilities for implementing 
information transfer along the QD-based bus, which exploit 
the global dynamics of the system. We shall devolve particu- 
lar attention to the constraints imposed on the scheme by the 
physical nature of the system itself. Later we shall discuss the 
robustness of the proposed schemes and compare the perfor- 
mances of our quantum bus against the use of a microcavity. 



n. QUANTUM COMPUTER HARDWARE DESIGN 

Self-assembled quantum dots are characterized by different 
properties, depending if we are considering an ensemble of 
in-plane QDs (the plane being orthogonal to the direction of 
growth) or an array of QDs stacked in the growth direction. In 
the following we will refer to the latter structure as 'Vertically 
Stacked Array' (VSA). 

VSAs have very distinctive characteristics. First of all 
they tend to form naturally due to strain propagation'. Sec- 
ondly self-assembled quantum dots are strongly confined in 
the growth direction, due to bandgap mismatch between the 
different semiconductors. This implies that (i) particle tun- 
neling between nearby, vertically stacked QDs can be more 
difficult, even for small QD-QD distances^ (ii) QDs can be 
very close together and (iii) interactions between vertically 
coupled quantum dots may be quite strong (see e.g. Ref|2|). 
The latter characteristic is very promising since quantum op- 
erations based on these interactions can be carried out in times 
much shorter than the relevant decoherence times^. Last but 
not least, due to the strain propagation, vertically stacked QDs 
tend to be noticeably different in size, which make them good 
candidates for energy selective addressing of specific excita- 
tions, a key ingredient in many QDs-based QIC proposals2iii. 

On the other side in-plane quantum dots can be produced 
in fairly uniform size by controlling the growth parameters: at 
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FIG. 1: Schematic sequence of vertical arrays (VSA) and in-plane 
arrays (IPA) of QDs. Different VSAs can be addressed by different 
laser spots. 



present samples can be already mono-disperse up to only few 
percent size variation. This uniformity will be most probably 
improved much further in the near future, especially since an- 
other important technological application, such as QD-based 
lasing, requires monodispersity'". In addition there is already 
a good control of the in-plane QD density (which can be var- 
ied between 10^ to 10^^ cm^^) and, at the moment, tech- 
niques, e.g. 'seeding', to produce more regular in-plane dis- 
tributions are proposed and testedii. 

The differences outlined above, between VSA and in-plane 
QDs, suggest us to diversify the use of the two categories, 
when thinking to the hardware structure of a quantum com- 
puter. The structure we envisage in this paper is composed 
by a series of VSAs - for example laying in the same plane - 
connected by chains of in-plane QDs, which we will refer to 
as In-Plane Arrays (IPAs, see Fig.fQ. 

This design allows for compact implementation of a quan- 
tum computer, since in principle it may confine the hardware 
in a single slab. This architecture satisfies the request of easier 
lateral access (e.g. for wiring) and manipulation of each QD 
forming the computer (e.g. for applying external fields to spe- 
cific dots). This design is also scalable, in the sense that the 
same basic structure (VSA plus related IPA) could be repeated 
at will, as well as the slab structure as a whole. 

We underline that the hardware structure we propose has 
yet another very advantageous feature. It can be built such that 
(i) each quantum register (i.e. any VSA) would be address- 
able by energy selective schemes (exploiting the natural diver- 
sity of stacked QDs), as required by many QIC schemes^'^'^, 
but, at the same time, (ii) different quantum registers (differ- 
ent VSAs) could be also discriminated spatially, e.g by near- 
field techniques. This can be achieved by making the quan- 
tum buses (i.e. the IPAs), longer or of the order of 1000 A( see 
Fig.G). In this way the different VSA can be spatially resolved 
by different laser pulses. This is a big improvement in respect 
to architectures which use energy selective addressing only: 
in our proposal even if the same laser color would in principle 
couple to one QD per VSA, we would be still able to address 
these resonant quantum dots separately. In general we could 
identify the QD sequence in each VSA with a short, well de- 
fined laser color sequence, and this sequence could then be re- 
peated (totally or partially) for different VSA. In principle the 
whole hardware structure could be built as a stacked sequence 
of mono-disperse QD planes, subsequently e.g. etched to form 
the VSA-IPA chosen structure. This combined spatial-plus- 



energy selective addressability strongly enhances the hard- 
ware flexibility and helps overcoming the problems related 
to schemes which use energy-selective addressing only (e.g. 
NMR-based schemes). 



III. QUANTUM-DOT QUANTUM-BUS POTENTIAL 
IMPLEMENTATIONS 

VSAs have been analyzed in details in many QIC 
proposals^'^'^; in this paper we will concentrate on possible 
ways of implementing a quantum-bus by using an IPA as hard- 
ware. 

We will discuss in details two potential implementations 
which exploit the global dynamics of excitons in the IPAs. In 
both cases we shall (i) consider that state-of-the-art in-plane 
QDs can be made mono-disperse up to few percents and (ii) 
discuss a way to fine-tune QD-QD interactions a posteriori. 



A. Implementation by global dynamics: general 
considerations 

Recently there have been interesting developments in the 
study of the global dynamics of a spin-chain'^-'-''''*''^ as a 
possible way of transferring quantum information across the 
chain. It is well known that the time necessary to trans- 
fer information along a spin or pseudo-spin chain in QIC 
sche mes bas ed on nearest-neighbor gating operations (e.g. see 
scales with the number of spins. There is though 
a widespread feeling that, when the transfer of information 
along a spin chain is based on the chain global dynamics, the 
transferring time can be independent from the chain length. 

We will explicitly show that, for real systems, the time tt for 
transferring a qubit of information across a spin chain still de- 
pends on the chain length (for our specific hardware the num- 
ber of QDs forming the IPA). In some cases though and for the 
same chain length, tt may be smaller when exploiting global 
dynamics than by using nearest-neighbor gating operations, 
and some other advantages may occur in the first case. 

To understand why this is so, we must take into consid- 
eration some very general constraints, which are due to the 
fact that any hardware systems is indeed a physical system. 
These constraints are often overlooked in more abstract de- 
scriptions of spin-chain global dynamics; when 'translating' 
the schemes into a physical system though, physical limita- 
tions must be taken into consideration. 

Let us consider a spin-chain global dynamics based on 
nearest-neighbor interactions Ji,i^i's which are modulated 
along the chain, 

J,,,+i = ./ • /(iV), (1) 

where J is a constant representing a typical energy scale of 
the system and f{N) is the function describing the coupling 
modulation along the chain. We assume the modulation to be 
a reasonably regular function and to be un-bounded, so that / 
will depend over the chain length N. 
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The physical constraints in applying this scheme to any 
concrete system are connected to the following observations: 
(i) in any physical system there will be an upper limit Jmax to 
the interaction strength, i.e. 



max{Ji^i+i} < J„ 



(2) 



and (ii) for practical purposes the transfer time U must be 
smaller than the typical decoherence times, i.e. 



U < Tde 



(3) 



We shall show that Eqs. (|2ji and (|3j (i) imply that the transfer 
time tt across the chain depends over the chain length N and 
that (ii) this sets a physical limitation to the maximum length 
of the chain we can consider. 

From Eq. (|2j and Eq. Q we have 



J-f{N) < J„ 



(4) 



Since J/h is the typical frequency of the system, we can as- 
sume that tt ^ {J/h)^^ . From Eq. @ we obtain 



f{N) 

J max 



< tt 



(5) 



which shows that (i) the transferring time has a lower bound 
and, more important, (ii) this lower bound (and hence the 
transferring time) depends on the chain length N. Using now 
Eq. (|3}, we obtain 



f{N) 

Jmax 



< Tde 



(6) 



which sets a physical upper limit to the mathematically un- 
bounded f{N) and hence to the chain length N. 

For chains in which the modulation of the couplings { J^} 
is bounded, the dependence of tt on the chain length is due 
at least to the constraint that information cannot travel faster 
than light. This is clearly a much weaker constraint, and, as 
we shall see below, our calculations show that, even for the 
extreme case of chains with un-modulated couplings, much 
stringent relations between tt and N usually occur. 

Among the various proposals related to spin chain global 
dynamics, we will first consider the one described in Ref. 12. 
This is a scheme for perfect transfer of information (a quan- 
tum state or two entangled states) across a spin chain, by ex- 
ploiting the global time evolution of the system. 

The requirement for perfect transfer is that the system is 
governed by a XY Hamiltonian, in which the non-diagonal 
couplings between nearby spins Ji,i+i must satisfy certain se- 
quences, e.g.'^ 



J^i{N-i), 



(7) 



where N is number of spins in the chain and i = 1, ...A^ — 1. 
Eq. Q represents an unbounded modulation of the coupling 
strengths of the type considered in Eq. Q, which is in addition 
symmetric in respect to the centre of the spin chain. The time 
requested for the perfect transfer isi^ 



(n+l)pfc 



(2n+ r 



flTI 



(8) 



n = 0, 1, 2..., and it is determined by the value of J. 

In particular let us consider the transfer time corresponding 
to n = 0, tipk = h'n/{AJ) and rewrite Eq. as 



J = J^ax/yJint{N/2) ■ [N - int{N/2)l 



(9) 



where we have optimistically assumed that at the center of the 
chain Ji^i+i = Jmax- Notice that the latter equation implies 
that 



t 



ipfe 



N h 

^ ^ m.a. 



(10) 



i.e. there is direct proportionality between minimum time of 
transfer and chain length N . If we consider a transfer across a 
qubit chain done by nearest-neighbor two-qubit gates, and the 
ideal situation in which all the physical couplings used in the 
gates are given by Jmax, we would have 



tlpk ^ N 



Jrr 



(11) 



Comparing the last two equations, we see that the actual gain 
in using the global dynamics of a spin chain, when the phys- 
ical condition of including a maximum strength for a cou- 
pling is taken into account, is a factor 2, coming from the 
symmetry of the modulation, but the transfer time still scales 
with the total number of spins N. Fig.|3] lower panel, shows 
tipk extracted from the simulation of spin chains with nearest- 
neighbor couplings following Eq. in respect to N. It 
confirms the direct proportionality between minimum trans- 
fer time and chain length. 

The concrete advantage of using global dynamics over two- 
qubit gating becomes then the possibility of simplifying the 
implementation, an issue that should not be underestimated, 
but may again depend over the specific physical character- 
istics of the system chosen as hardware. At the very least, 
exploiting global dynamics lifts the necessity of driving every 
single two qubit operation and this consequently avoids poten- 
tial sources of errors; but what is really important when con- 
sidering a spin chain global dynamics as a way to implement 
a quantum bus, is to maximize the simplicity of the device. 

The fact that, as shown above, for physical realizations, N 
cannot be arbitrarily long, implies that the properties relevant 
to many practical implementations will be the ones of medium 
or relatively short chains, while previous literature has mainly 
focused on the large-N limisi^*i^. 

With this in mind, in the following we shall focus on chains 

characterized by iV ^ 10, both underlining some specific gen- 
eral properties and referring to a specific physical system - the 
IPA-chain formed by QDs with excitons as pseudo-spins. 



B. Forster effect 

We will consider implementations of quantum buses in 
which the coupling between nearby in-plane QDs is provided 
by the Forster effect'^ '^. Due to this coupling, the excitation 
corresponding to a correlated electron-hole pair (exciton) can 
be transferred from one QD to its neighbor. 
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This effect has been observed as a powerful motor for exci- 
tonic transfer in colloidal self-assembled quantum dots'^. The 
measured strength of the Forster coupling in CdSe QDs is of 
about 0.3 meV, but it is believed that in structurally optimized 
systems the interdot energy transfer can approach picosecond 
time scales, i.e. meV coupling strength'^. Being the structures 
of interest to this paper (self-assembled QD) of similar size in 
respect to colloidal QD and being typical values for QD-QD 
interactions of few meV^, we will assume in our calculations 
1 meV as a reasonable upper limit for Forster interaction in 
self-assembled QDs. 

In the envelop function approximation, the Forster coupling 
between an exciton in QD^ and an exciton in QDj, Vpij, is 
given by— 

- l^(<i-«>-Ru)(<r, >-Ry)] (12) 

where 

= 1 rl:l{v)rh{v)d\ (13) 

is the overlap between the hole and electron envelop functions 
in QDi, Ri j is the distance between QD, and QDj, and e < 
Yi > is the dipole moment due to the atomic part of the wave 
function in QD;. 

We underline that no particle tunneling is involved in the 
Forster process, but only energy transfer: it depends in fact 
on Oi, the overlap of electron and hole wave-functions in the 
same QD. In tunneling, the relevant quantity is instead the 
overlap of wave-functions belonging to different QDs. Due to 
energy conservation, this coupling is most effective between 
very similar QDs. When QD chains are involved though, we 
can picture a dynamical Forster-based process between far- 
away but similar QDs positioned along the chain. The simi- 
larity of the first and last QD in the chain ensures that the total 
process of transferring the excitation across the IPA remains 
energy preserving. 

Labeling with |0)i (| l}i) the absence (the presence) of a sin- 
gle ground-state exciton in QDj, the Hamiltonian for a chain 
of N QDs, in which we consider at most single exciton occu- 
pancy of the IPA, is given by 

N 

i 

+ J2 {yF^,^+l\^)^{0\^ ® |0),+l(ll,+ l + h.C.) ,(14) 

where we have assumed no single particle tunneling. Contri- 
butions behind nearest-neighbors are second order, since the 
Forster effect scale as i?? , and have been neglected. 



C. Implementation by global dynamics using a QD chain: first 
method (modulated coupling) 

For concreteness, in the following we shall consider as 
chain an IPA, and as (pseudo) spins degrees of freedom the 
presence or absence of a ground state exciton in each QD-. 
The transfer of quantum information across the "spin chain", 
will then correspond to the transfer of an excitonic state across 
the IPA, by means of exploiting the global excitonic dynamics 
of the system. The general conclusions on the IPA can though 
be easily transferred to other types of physical hardware. 

The method we are going to discuss in details represents 
one of the possible schemes for transferring information be- 
tween separate VSAs in the computer architecture previously 
outlined. To be able to address different VSAs with inde- 
pendent laser pulses, we require the length of the IPA to be 
longer than 1000 A. A reasonable distance between the center 
of nearby QDs is at least of the order of 150 — 20GA. This 
implies that the IPA must be formed by at least 5-7 QDs. The 
question is now if excitonic transfer through such a chain can 
be performed in a reasonable time and with reasonable accu- 
racy. 

In our system the (pseudo) spin degrees of freedom are rep- 
resented by the presence (|l)i) or absence (|0)i) of a ground 
state exciton in each QD^, where i — 1,N labels sequentially 
the QDs in the chain. This is the same encoding proposed in 
Ref.H The Hamiltonian of the system is given by Eq. il4\ . 

We shall refer to the scheme presented in Ref. 12 as the 
Modulated Coupling (MC) scheme. In respect to the XY 
Hamiltonian required by Ref. J2ai Eq. ( I14t presents diagonal, 
in general not-equal, terms Ei and XY non-diagonal terms, 
provided by the Forster coupling Vpij, which in general do 
not satisfy Eq. 0. To obtain a good approximation to the 
'mathematical' Hamiltonian, we need at least to engineer the 
non-diagonal terms into the required sequence Eq. Q. From 
Eq. il2\ . we see that we can tune the value of the Forster cou- 
pling by (i) engineering the value of the interdot distance Rij 
or (ii) by tuning the value of the overlap integrals Oi. Since 
roughly Vpij ^ l/R^y to engineer Vpij according to Eq. 
over an IPA of 5-7 QDs with Ri - 150 - 200A, implies to 
be able to experimentally control Rij up to few Amstrongs, 
which is at the moment experimentally very challenging. A 
more flexible approach, is to vary, a posteriori, the electron- 
hole function overlap Oi. This is doable by applying a suitable 
(but static) in-plane electric field to each QD, as described in 
Ref. 2. The field is to be applied in the in-plane direction or- 
thogonal to the IPA's direction (see Fig. |5} . The advantage 
of this scheme is to allow also for corrections of possible er- 
rors occurred in the growth process, since the applied electric 
field can be varied a posteriori, after the growth process is 
completed. 

Let us consider a 'uniform' initial condition, in which 
Vpi.i+i = Vfjj+i for every i, j. If TV = 5-7and {VFj,i+i} 
has to satisfy Eq. Q, the maximum correction required over 
the coupling strength Vpij is of about 30%. This corresponds 
to a 20% variation over the overlap integrals O^'s. We need 
now to estimate the required electric field, in order to asses the 
feasibility of the scheme. If we assume Gaussian electron and 
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FIG. 2: Sketch of an IPA in which an external static electric field T\ 



has been applied to every QDi 



hole wave-functions, the relation between Oi and the distance 
d between the wave-function centers is 




1 



exp 



9 X2 I \2 



(15) 

where \te{h) ^"e the characteristic length of the wave- 
functions in the t — x,y,z directions. If we assume similar 
dots (a reasonable request for an IPA), we obtain 



'2(AL + AL) 



In 



O 



np 



(16) 



where Onp, is the unperturbed value of the overlap integral, 
i.e. the value it assumes for the central(s) QD(s) which is 
coupled to the neighbors by the coupling constant with the the 
highest value (see Eq. (0). On the other side di can be also 
related to the static electric field Ti applied to QD^ by^ 



eTiiil/nieUjl^) + {l/mtujl^,)]. 



(17) 



Assuming typical GaAs parameters (irih — 0.34mo, nie = 
0.067mo, hujh ^ 20 meV, hoje ~ 40 meV), and using 
Eqs. il6\ and illi . we see that a variation of 20% over the 
overlap integral Oi can be produced by a field of the order of 
30KeV/cm, perfectly reasonable from the experimental point 
of view. We can then assert that the scheme we are going to 
describe more in details could in principle be implemented by 
state of the art technology. 



D. Modulated coupling: simulations and errors 

In this section we shall discuss and simulate the excitonic 
dynamics in IPAs and discuss the effect of possible sources 
of errors for the transmission of the excitonic state across the 
IPA. Because of the generality of most of the results, we will 
indicate the couplings as Ji.i+i, keeping in mind that when 
the considered hardware is an IPA the coupling is given by the 
Forster effect Eq. J12t . We will consider in our calculations 
realistic QDs parameters and in particular GaAs-based IPAs 
with < 9 and Jmax < 2meV, Jmax the upper limit for the 
coupling between nearby QDs. 

Let us first of all consider the case in which the condi- 
tions in Ref. are perfectly satisfied. Fig. |3l (upper panel) 
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FIG. 3: Upper panel: {10000|*(t)> (dashed-dot line) and 
(000011\E'(f)) (solid line) vs time for N=5. 'Ipk' marks the first 
resonance position. Lower panel: fip^ vs N. 



shows a simulation of the dynamics of the system for Jmax = 
ImeV and an IPA of 5 QDs. The state at time i = is 
[■(/"(O)) = jlOOOO), i.e. single exciton in QDi. The projec- 
tion of the state of the system at time t on the initial 
and the desired final state (i.e. exciton in QD^r, the last QD), 
(10...0|i/'(0) andjO... 01 !?/;(*)) respectively, are plotted. As 
predicted in Ref.UZ, there is periodic perfect transfer between 
the two states, the first transfer occurring at tipk ~ 1.3ps. 
The lower panel presents the plot of the numerical value ob- 
tained for tipk when the length of the IPA is changed. As 
can be seen, the linear dependence of Eq. jlO> is recovered. 
Notice that, even for = 9, the transferring time remains 
very short, tipk ~ 2.3ps, confirming that, even if there is di- 
rect proportionality between N and tipk, at least for relatively 
short chains, the effect of this proportionality is not dramatic. 

Let us compare these results with the use of microcavities, 
which is the alternative route for transferring information at 
'long distance' in QD-based quantum computers. If we con- 
sider a comparable transfer done by mean of a microcavity, 
i.e. with a coupling of 1 meV between QDs and cavity, the 
minimum time necessary for the transfer is greater than 100 
ps^. Extrapolating from Fig. |3j lower panel, and assuming 
the relation Eq. ( I10> . a transferring time tipk — lOOps would 
correspond to a transfer over an IPA with N « 400, i.e., as- 



suming Rij — 150 ~ 200A, over a distance of 6 S/xm. We 
see than that, from the transfer-time point of view, the use of 
microcavities becomes advantageous only when really long 
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distances, of the order of several microns, are involved. For 
these distances the design and modulation of long spin-chains 
(long IPAs) may also become too complicated. 

We notice also that a phase shift is originated at the end of 
the transfer, depending on the number of QDs N—. In par- 
ticular when — 1 is an exact multiple of 4, no phase shift 
is produced-2. Since, once the length of the chain is known, 
this fixed shift can be corrected by a simple single qubit phase 
rotation operation done after the transfer, we will neglect the 
effects of this shift in the following discussion. 

Let us focus now on our specific system, the IPA. If we 
consider as a starting point a chain of similar QDs regularly 
spaced, and approximate the electron and hole wave-functions 
as Gaussian, we can estimate the Stark shift of the excitonic 
level of QDi in respect to the central, unperturbed, QD„p pro- 
duced by the application of the fields {JT^ } to satisfy Eq. Q 
(see Fig.|2ji. This is given by 



1 



\E,, 



■In 



(18) 



where Jmax = max{Ji] is the value of the coupling between 
QD„p and its neighbors. The ratio Ji/ Jmax can be calculated 
from Eq. Q. For a given semiconductor, Eq. ( I18> depends 
on the ratio LOex.i/^hx.i, but for the relatively short IPAs we 
are considering, we see that, even when u)ex.il^hx,i ~ 3 and 
huiex.i = 20 — 40 meV, \Ei\ > ImeV. When we consider the 
approximate sequence {Ei} as given by Eq. ( I18> and let the 
system evolve according to the Hamiltonian in Eq. J14> . our 
simulations show that the excitonic transfer from the first to 
the last QD may be achieved, but depending on the param- 
eters chosen. In Fig. 0] we compare the case of N=5 (upper 
panel) and N=7 for the same structural parameters (see fig- 
ure caption) and |V'(0)} = |10...0). Even when the transfer 
is achieved (N=5), the periodic pattern becomes more com- 
plex and other excitonic states of the system are populated 
during the transfer itself (e.g. the probability of occupying 
state jOOlOO) is shown). For N=7, over the same interval of 
time, no relevant transfer across the chain is achieved, while 
the exciton remains mainly localized in QDi. 

We can think of improving this result for relatively short 
IPA by requiring careful tuning of at least part of the distances 
Rij between QDs and/or of other QDs parameters during the 
growth process, in order to obtain as a starting point a se- 
quence of Ji which matches more closely the requirement of 
Eq. 0. This implies the application of smaller static elec- 
tric fields and consequently smaller Stark shifts. The energy 
shifts will now be in general randomly distributed, depending 
on how well the single coupling Ji^i+i satisfies Eq. Q. 

We have tested the robustness of this scenario by calculat- 
ing the fidelity of the transfer across a QD chain in which 
the sequence {Ei} is produced randomly, with the constraint 
Ei < SE. The fidelity of the transfer for each set {Ei} is 
calculated according to Ref.ll3las 



{F)bi = 



|/Ar|cOS7 I/atP 



1 



(19) 
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with /jv = (0....1|exp(-ii?i)|1....0), H given by Eq. (EJ 
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FIG. 4: Upper panel: |(10000|i/'(t))p (solid line), | (00001|i/)(t)> ^ 
(dashed-dot line) and | (00100|i/'(f)> ^ (dotted line) vs time, for 
Jmax = 1 meV, hujex.i = 20 meV and ujex.i/i^hx.i = 3 and GaAs 
parameters. 

Lower panel: |(1000000ii/'(i))i^ (solid line), \{000000l\ip{t))f 
(dashed-dot line) and j (0001000|^/)(t)) ^ (dotted line) vs time. Same 
parameters as upper panel. 



and 7 = arg}/^}. In the following we will set cos 7 = 1 
since we can always correct the phase shift after the transfer 
has occurred, as mentioned before. Notice that the fidelity in 
Eq. ( II9I 1 includes the averaging of the initial state of QDi over 
the whole Bloch sphere'^. The obtained result will be then 
averaged (when appropriate) over different random configu- 
rations {Ei} ({F)Bi,cf)- In this respect, when 5E ^ Jmax it 
is important to perform this average over some thousand con- 
figurations, in order to obtain reliable results. 

Results for IPAs with N=5 and N=7 are plotted in Fig. [S] 
{F) Bi,cf is plotted in respect to the value of Jmax, for two dif- 
ferent values of SE. The behavior for the two values is similar, 
the relevant parameter being the ratio SE/Jmax In particular, 
when SE « Jmax, the fidelity is still greater than 90%, but it 
decreases drastically for increasing SE/Jmax- The variation 
with the chain length is instead not substantial, though shorter 
chains perform better Notice that in Ref.Ts'the opposite trend 
is predicted, i.e. longer chains more robust than shorter ones. 
This is due to the fact that in their simulations they use dif- 
ferent constraints: they vary N for fixed time transfer tipk, 
i.e. for fixed J. In other words the maximum value for the 
coupling between nearby spins is not fixed, and Jmax is let 
increasing with the number of spins N. 
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FIG. 5: Average fidelity {F)Bi,cf vs Jmax, for 5E = 0.1 meV and 
SE — 1 meV and chain length of 5 and 7 QDs, as labeled. The 
vertical dashed lines indicates the values for which Jmax ~ SE. 



When SE ^ Jmax, the fidelity of the transfer (F) bi from 
first to last QD largely varies according to the specific set 
{Ei}. In addition the presence of 0} modifies the 

time at which (F) bi is maximum (the same occurs for the 
second method we shall consider, see inset of Fig.|8|i. In this 
respect a possible way to improve the fidelity would be to ini- 
tially measure the optimal transfer time for the specific set 
{Ei} corresponding to the sample and then tune the time for 
the transfer operation accordingly. For IPAs the main contri- 
bution to {Ei] is given by growth and/or applied static elec- 
tric fields, so that measuring the optimal transfer time for the 
specific hardware piece might recover a significant part of the 
fidelity. Nowadays tests on single electronic components are 
done routinely to check their specifications; we foresee that 
analogous testing procedures will be needed for the various 
hardware components of a quantum computer. 



E. Implementation by global dynamics using a QD chain: 
second method (flat coupUng) 



We have stressed in Sec. IIII Al that the main advantage of 
using global dynamics should be to simplify the implementa- 
tion of the computation. For implementing the MC method 
with IPAs, carefuU tuning of the QD properties is required 
both in the growth phase and a posteriori, with the risk that 
large values of Ei may result, affecting the fidelity. In order 
to simplify the implementation and overcome some possible 
sources of error, we will consider in this section the global 
dynamics of a spin chain, when all couplings between nearby 
spins are the same, i.e. Jk = Ji, and Ek = Ei for every 
i, k. We shall refer to this as "flat coupling" (FC) scheme. As 
underlined in Sec. HI] these requirements should be easier to 
match, at least considering implementations which use IPAs 
as hardware, since mastering these QD properties is already 
on the agenda of many experimental groups"^". 

We want to address the questions (i) if transfer of a quantum 
state across the chain by exploiting global dynamics is still 



possible and with which degree of accuracy and (ii) how the 
results are affected by imperfections. 

We consider again relatively short spin chains, with N — 
1. ...9, and simulate the evolution of the system according to 
the Hamiltonian Eq. J14t . A similar problem have been an- 
alyzed in Ref. |3 but under different constraints and mainly 
considering the limit of long spin chains, so that properties 
of the relatively short chains which are relevant to the present 
proposal have not been discussed. 

The inset of Fig. |6l shows the time evolution of a N=5 chain 
with Ek = 0, k — 1,N and Jmax = ImeV. The upper 
panel presents the fidelity (F) bi in respect to time; the lower 
panel shows the real (solid line) and the imaginary (dashed- 
dot line) components of the projection (00001 in re- 

spect to time, with initial conditions |*(0)) = |10000). As 
can be seen, the time evolution as a whole presents a periodic 
behavior (within the numerical error), centered on the peak la- 
beled with 'C . The period is proportional to 1/ Jmax- Within 
this period many high-fidelity resonances are present, some of 
which correspond to a fidelity {F)bi ^ 1 (e.g.(_F)^j = 0.981, 
= 0.996 and = 1.00). The behavior just dis- 

cussed represents a typical pattern for all the spin-chains here 
analyzed. In particular, for < 9, the very first resonance 
(labeled 'A' in the inset of Fig.|6jl corresponds always to a fi- 
delity {F)bi > 94% (Fig.0 main panel). In general chains 
with an odd number of QDs N have an overall period shorter 
than the corresponding (N-1) QD chains. Interestingly the 
phase of the projection (00001 1 'J (t)) follows for FC the same 
dependence on N as for MC, i.e. it is given by (— i)^^^. 

In the inset of Fig.0 the time <^ at which the resonance 'A' 
occurs, is shown in respect to N. Both and the correspond- 
ing fidelity {F)^i present a weak linear dependence with the 
length of the chain N, with though still below 2 ps for N=9 
(inset of Fig.Q. We stress that for the applications we are in- 
terested in, we need a finite, relatively small number of QDs, 
so these results are very encouraging. 

In all cases analyzed, at least one of the resonances follow- 
ing peak 'A' corresponds to an even higher fidelity (e.g. B 
and C in Fig.|6j for N=5). The possibility of accessing these 
stronger resonances as a mean of a better information trans- 
fer along the chain, will depend on the values of Jmax, and 
of the typical decoherence times of the system, since, as al- 
ready mentioned, the times corresponding to the different res- 
onances scale as 1/ Jmax- 

We shall now analyze how the different resonances are in- 
fluenced by imperfections. As for MC, we have considered 
sets {Ei} with Ei\ randomly distributed and Ei < 5E^'i. 
An average over an appropriate number of such configura- 
tions is taken to calculate the fidelity {F)Bi,cf- We present 
the results in Fig. |5] where the average fidelities {F)gi ^y, 
K = A, B, C, calculated at the times Ia, tB and tc at which 
the resonances occur, respectively, are plotted against SE. Our 
results show that the resonances B and C are less 'robust' in 
respect to the first resonance (A), i.e. they tolerate a lesser 
amount of imperfections in the flat distribution of the energies 
Ek- In fact for very small SE/ Jmax, {F)^'i^cf * while 
{F)^i = 0.98. {F)g'flj: though decrease rapidly already 
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FIG. 6: {F)Bi,cf vs SE for three different resonances (labeled A,B 
and C), for N=5 and Jmax = 1 meV. The dashed vertical line corre- 
sponds to the value of Jmax- Inset: upper panel: {F) bi as a function 
of time for SE — 0, N=5 and Jmax = 1 meV. Lower panel: real 
(solid line) and imaginary (dashed-dot line) part of (00001|^(t)) for 
1^(0)) = 1 10000) and same parameters as upper panel. 
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FIG. 8: {F) Bi,cf vs SE for MC scheme at tipk (solid line) and FC 
scheme at tA (dashed-dot line). N=5, Jmax ~ 1 meV. Inset: (F) bi 
vs time for FC scheme and same parameter as main panel. The x-axis 
upper limit corresponds to tA- 



when using IPAs as hardware), it is convenient to opt for the 
simpler FC scheme. The inset of Fig.|S]presents (F) Bi{t) for 
FC in respect to time, for t < Ia and 5E — 1 meV. Each 
curve corresponds to a different random set {Ei}. The plot 
shows that {F)^i is spread over a large interval which in- 
creases with increasing 6E (not shown). Two things should be 
noticed, the first is that the fidelity corresponding to a specific 
{Ek] configuration can be (clearly) greater than the average 
fidelity {F)Bi,cf, the second that the maximum the fidelity 
for the transfer may occur at times different from tA- Similar 
conclusions are reached when analyzing (F) bi {t) in the MC 
scheme. These considerations imply that (i) operatively it will 
be important to carefully measure the sample characteristics in 
order to optimize the transfer time, and/or, as possible alterna- 
tive, (ii) specific sets {Ei} could be designed to optimize the 
fidelity. 



FIG. 7: {F)bi vs N. Inset: tA vs N. Parameters as inset of Fig.|S| 

for SE/Jmax ~ 0.1, while {F)^^ remains basically unper- 
turbed. For 5E/J^ax = 1, {F)Bfcf 

drop to about 0.7, while 
{F)gi ^ 0.94, still a very significant fideUty. 

F. Comparison between modulated and flat coupling 

We shall now compare MC and FC schemes. In Fig.|8]we 
compare {F)"^i (dashed-dot line), with the average fidelity 
{F)Bi.cf calculated at tipk (see Eq. (|8|i) with J^'s given by 
Eq. Q (solid line). They are plotted in respect to SE, which 
is varied up to 2 meV, for N=5 and Jmax = 1 meV. We notice 
that for SE/ Jmax ~ 1-3, the fidelity of MC method becomes 
smaller, confirming that, if SE ^ Jmax is to be expected (e.g. 



FC scheme can be seen as representing a very bad imple- 
mentation of the MC scheme, in which every trace of the se- 
quence Eq. has been washed out. In this respect the oc- 
cuiTence of the very strong first resonance 'A' - resonance 
which coincides, time and strength, with the first transfer of 
the MC scheme when iV < 3 - suggests that, for relatively 
short chains, the MC scheme is very robust in respect to im- 
perfections in the J^'s sequence Eq. 0. This robustness de- 
creases with N, as the strength of {F)gi does (see Fig.El^"- 
This interpretation also helps to understand why (i)the same 
phase dependence on N is observed in the FC and in the MC 
scheme for the projection (00001|^(i)) and (n)tA has a lin- 
ear dependence on N similar to tipk (compare Fig.|3] lower 
panel, and inset of Fig.0. 

Our results suggest, for simplicity of implementation and 
especially if it is possible to have access to the resonances 
corresponding to the highest fidelities (B and C in Fig.|6jl the 
FC as the best scheme, at least when using IPAs as hardware. 
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IV. SUMMARY AND CONCLUSIONS 

We have proposed a potential implementation for a quan- 
tum computer where both the computation and the commu- 
nication between far away qubits utilize exclusively quantum 
dot arrays. We have discussed how the different characteris- 
tics of in-plane and stacked arrays of QDs can be exploited 
to this purpose, the first as quantum buses and the second as 
quantum registers. The structure we propose is scalable and 
avoids the use of microcavities. 

In the second part of the paper we have focused on the prop- 
erties of the global dynamics of 'relatively short' spin chains, 
such chains been the most relevant in respect to the physi- 
cal problem at hand. By mapping ground state excitons into 
(pseudo) spins, this dynamics can be used as a way to im- 
plement long distance communication between far away QDs 
connected by an in-plane array of QDs. We have discussed 
and compared two different potential schemes, the first based 
on the coupling between nearest-neighbor spins modulated 
along the chaini^, the second based on a flat, uniform cou- 
pling. In particular, in the case of flat coupling, we have shown 
that the system presents a behavior overall periodic in time, 
containing a series of resonances corresponding to a fidelity 
of the state transfer across the chain very close to 100%. The 
very first of these resonances corresponds always to a high fi- 



delity, larger than 94% for N < 9. This fidelity presents a 
linear dependence on the spin chain length N, but such a de- 
pendence is relatively weak. According to the system physical 
characteristics (maximum value of the nearest-neighbor cou- 
pling and decoherence times) one or all of these resonances 
will be accessible experimentaUy and exploitable as a mech- 
anism for transferring information. The effects of imperfec- 
tions on these resonances has also been analyzed, showing 
that the very first one is very robust. 

It is clearly possible to design many more schemes to 
implement information transfer along in-plane QD chains - 
schemes either based on the global dynamics of the chain or 
on nearest-neighbor two-qubit gating. A systematic analysis 
of all possible schemes goes beyond the scope (and the possi- 
bilities) of this paper. We hope that the two potential imple- 
mentations here discussed represent a good starting point for 
a stimulating discussion on the all-QD-based hardware struc- 
ture we have proposed. 

Acknowledgments 

We ai-e grateful to M. Feng, N. Datta, M. Christandl, A. 
Bychkov, TP. Spiller, and B.W. Lowett for stimulating and 
fruitful discussions. 



* Electronic address: |ida500@york.ac.uk| 

' Pierre M. Petroff, Axel Lorke, and Atac Imamoglu, Physics Today 

54, issue 5, 46 (2001). 
^ E. Biolatti, I. D'Amico, P. Zanardi and F.Rossi, Phys. Rev. B 65, 

075306 (2002). 

' S. De Rinaldis, I. D'Amico, E. Biolatti, R. Rinaldi, R. Cingolani, 

and R Rossi, Phys. Rev. B 65, 081309 (2002). 
" A. Nazir, B.W. Lovett, S.D. Barret, TP Spiller and G. A.D. Briggs, 

Phys. Rev. Lett. 93, 150502 (2004); A. Nazir, B.W. Lovett, and 

G.A.D. Briggs, Phys. Rev. A 70, 052301 (2004) 
^ M. Feng, 1. D'Amico, P. Zanardi and F. Rossi, Phys. Rev. A 67, 

014306 (2003). 

* M. Feng, I. D'Amico, P. Zanardi and F. Rossi, Europhys. Lett. 
66,14(2004). 

' E. Pazy, E. Biolatti, T. Calarco, I. D'Amico, P. Zanardi, F. Rossi, 
R Zoller, Europhys. Lett. 62, 175 (2003). 

* M.S. Sherwin, A. Imamoglu, and T. Montroy, Phys. Rev. A 60, 
3508 (1999) 

' S. De Rinaldis, I. D'Amico, and F. Rossi, Phys. Rev. B Phys. Rev. 
B 69, 235316 (2004). 

A. Moritz,R. Wirth, A. Hangleiter, A. Kurtenbach, and K. Eberl, 
Appl. Phys. Lett., 69, 212 (1996); Y K Su, S J Chang, L W Ji, C 
S Chang, L W Wu, W C Lai, T H Fang and K T Lam, Semicond. 
Sci. Technol. 19 No 3, 389 (2004);PR. Edwards, R.W. Martin, 
I.M. Watson, C. Liu, R.A. Taylor, J.H. Rice, J.H. Na, I.W. Robin- 
son and J.D. Smith, Appl. Phys. Lett., 85, 4281 (2004). Gao, M. 
Buda, H. H. Tan, C. Jagadish, Electrochem. Solid-State Lett. 8, 57 
(2005) 

G. Medeiros-Ribeiro, R. L. Maltez, A. A. Bemussi, D. Ugarte, 



and W. de Carvalho, Jr., Journal of Appl. Phys. 89, Issue 11, 6548 

(2001) ; W. Ye, S. Hanson, M. Reason, X. Weng, and R.S. Gold- 
man, J. Vac. Sci. Technol. B 23, 1736 (2005); X. Weng, W. Ye, 
S.J. Clarke, R.S. Goldman, V. Rotberg, A. Daniel, and R. Clarke, 
J. Appl. Phys. 97, 064301 (2005); Z. M. Wang, K. Holmes, Yu. 
I. Mazur, and G. J. Salamo, Appl. Phys. Lett. 84, Issue 11, 1931 
(2004). 

M. Christandl, N. Datta, T. C. Dorlas, A. Ekert, A. Kay and A.J. 

Landahl, Phys. Rev. A 71, 032312 (2005) 

S. Bose, Phys. Rev. Lett. 91, 207901 (2003) 

A. Woicik, T. Luczak, P Kurzynsky, A. Grudka, T. Gdala and M. 

Bednarska, quant-ph/0505097 , L.F. Santos quant-ph/0505072 

G. De Chiara, D. Rossini, S. Montangero and R. Fazio, 

quant-ph/0502148 

R.S. Knox and H. Van Amerongen, J. Phys. Chem. B 2002, 5289 

(2002) . 

" S.A. Crooker, J.A.Hollingsworth, S. Tretiak, and V.I. Klimov, 
Phys. Rev. Lett., 89, 186802 (2002). 

'^^ B.W. Lovett, J.H. Reina, A. Nazir, B. Kothari and G.A.D. Briggs, 
Phys. Lett. A, 315, 136 (2003). 

" By applying external electric fields, imperfections on the coupling 
constants may be corrected and partially transformed into imper- 
fections on the {Ei} set (as seen in the first part of the previous 
section), so for our purposes it is sufficient to study the effects of 
imperfections over the energy set. 

This conclusion is different from the one in Ref. 15 for the same 
reasons as for the effect of imperfections over the set {Ei}, 



